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FAST APPROACH TO THE TRACY- WIDOM LAW AT 
THE EDGE OF GOE AND GUE 

By Iain M. Johnstone* and Zongming Ma* 
Stanford University and University of Pennsylvania 

, We study the rate of convergence for the largest eigenvalue dis- 

■ tributions in the Gaussian unitary and orthogonal ensembles to their 

Tracy- Widom limits. 

^ We show that one can achieve an 0{N~^^^) rate with particular 

choices of the centering and scaling constants. The arguments here 
also shed light on more complicated cases of Laguerre and Jacobi 
ensembles, in both unitary and orthogonal versions. 

Numerical work shows that the suggested constants yield reason- 
' able approximations even for suprisingly small values of A'^. 

'. 1. Introduction. The celebrated papers of Tracy & Widom (1994, 1996) 

^ ' described the hmiting distributions of the largest eigenvalues of the Gaus- 

sian unitary and orthogonal ensembles (GUE and GOE) respectively. The 
purpose of this article is to show that an appropriate choice of centering 
and scaling allows us to establish a rate of convergence in these results, and 
further that this rate can be understood as "second order", being 0{N~'^/'^) 

■ rather than the 0{N~^/^) that would otherwise apply. 
The Gaussian ensembles refer, as is usual, to eigenvalue densities of x = 

(xi, . . .,xn) given by 

N 

f{x) = CNpJ{e-P^"'''J{\xi-Xj\P, 

i=l i<j 

with /3 = 1 corresponding to GOEtv and /3 = 2 to GUEtv, the subscript 
being shown only when clarity dictates. The corresponding matrix models 
' specify that / is the density of the eigenvalues x of a symmetric or Hermitian 

random matrix M with independent entries on and above the diagonal, 
whose density function is given by 

5(M) =4^exp{-(/3/2)trM2}. 
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Our principal rate of convergence results follow. The Tracy- Widom dis- 
tributions are denoted Fp{s) for /3 = 1,2. 

Theorem 1. Let denote the largest eigenvalue of a sample from 
GUEn, and 

fiN = V2N, TN = 
Given sq, there exists C = C{sq) such that for s > sq, 

- fiN)/rN <s}- F2{s)\ < CiV-2/3e-^ 

Theorem 2. Let denote the largest eigenvalue of a sample from 
GOEjsf^i, with + 1 even, and 

(1) m = V2N + 1, TN = 2"i/2iV-V6. 
Given sq, there exists C = C{sq) such that for s > sq, 

(2) - piN)/TM < s} - < CN-^'^e~'l\ 

We use index + 1 (rather than A^) because of a key formula relating the 
Gaussian orthogonal ensemble GOEat+i to the Gaussian unitary ensemble 
GUEjv, (36) below. The centering and scaling constants carry subscripts N 
rather than + 1 for this reason. 

Our interest in these results is threefold. First, they provide the simplest 
case of a class of such 0{N~'^/^) convergence results for the classical orthog- 
onal polynomial ensembles — the other two being the Laguerre and Jacobi 
ensembles — in both orthogonal and unitary versions. These results are of 
interest in statistics because they show that the Tracy- Widom approxima- 
tion is accurate enough to replace exact evaluation of the finite LOE and 
JOE probabilities for many applied purposes where highly accurate values 
are not necessary (Johnstone 2009). The results of this paper focus on the 
corresponding phenomenon for the simplest case of GUE and GOE. Since 
the LOE and JOE proofs are lengthy analyses with Laguerre and Jacobi 
polynomial asymptotics respectively, this paper outlines the approach in 
the simplest case. 

Second, our interest was stimulated by Choup (2009), which provided the 
leading terms in an Edgeworth expansion of the largest eigenvalue distribu- 
tion of GOE, and remarked that the N~^/'^ correction term does not vanish 
in GOE. As our earlier results on 0{N~'^/^) convergence for LOE and JOE 
would suggest a similar 0(A^~^/^) property for GOE with a suitable specific 
centering, it seemed, therefore, of interest to verify the conjecture in this 
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setting. Although we subsequently learned of an error in the argument of 
Choup (2009) (private communication), it was an important stimulus for 
this work. 

Third, we find it of interest that adjustment of and tn to secure 

0(iV-2/3) 

convergence yields an approximation that is adequate — for some 
purposes — for suprisingly small values of N . 

To illustrate, first in GUE, Table 1 shows the exact probabilities < 
^J'N + ''"7VS2o} for quantiles S2a of the limiting F2 distribution, computed 
using the finite GUE function provided in the Matlab toolbox RMTFredholm 
(Bornemann 2011). 



N 


MAT 


0.01 


0.05 


0.1 


0.3 


0.5 


0.7 


0.9 


0.95 


0.99 


2 


2.000 


0.026 


0.087 


0.149 


0.359 


0.549 


0.732 


0.913 


0.958 


0.992 


5 


3.162 


0.017 


0.068 


0.125 


0.331 


0.526 


0.717 


0.907 


0.954 


0.991 


10 


4.472 


0.014 


0.061 


0.115 


0.319 


0.516 


0.711 


0.905 


0.953 


0.991 


25 


7.071 


0.012 


0.056 


0.108 


0.310 


0.509 


0.706 


0.902 


0.951 


0.991 


50 


10.000 


0.011 


0.054 


0.105 


0.307 


0.506 


0.704 


0.902 


0.951 


0.990 


75 


12.247 


0.011 


0.053 


0.104 


0.305 


0.504 


0.703 


0.901 


0.951 


0.990 


100 


14.142 


0.011 


0.052 


0.103 


0.304 


0.504 


0.702 


0.901 


0.951 


0.990 


200 


20.000 


0.011 


0.051 


0.102 


0.303 


0.502 


0.701 


0.901 


0.950 


0.990 


500 


31.623 


0.010 


0.051 


0.101 


0.301 


0.501 


0.701 


0.900 


0.950 


0.990 



Table 1 

GUE approximation. For each percentile a shown in the top row, let S2a be the quantile 
F2{s2a) = a, and sjvq ~ ^jv + tnS2ci for /xjv and tn specified in Theorem 1. Table 
entries are -P{a;(i) < snc} computed using Bornemann's code for E2"\0; [sjva,oo)). 



In fact, our proof suggests a slightly different centering value, = 
{■\/2N — 1 + \/2N + 1) /2, which differs from \/2N in relative terms by only 
0{N~'^). However, Table 2 shows an observable improvement at very small 
values of N. 



N 




0.01 


0.05 


0.1 


0.3 


0.5 


0.7 


0.9 


0.95 


0.99 


2 


1.984 


0.025 


0.083 


0.143 


0.349 


0.538 


0.723 


0.909 


0.955 


0.992 


5 


3.158 


0.017 


0.067 


0.123 


0.328 


0.523 


0.715 


0.906 


0.954 


0.991 


10 


4.471 


0.014 


0.061 


0.115 


0.318 


0.515 


0.710 


0.904 


0.952 


0.991 


25 


7.071 


0.012 


0.056 


0.108 


0.310 


0.508 


0.706 


0.902 


0.951 


0.990 


50 


10.000 


0.011 


0.054 


0.105 


0.306 


0.505 


0.704 


0.901 


0.951 


0.990 



Table 2 

GUE approximation: as for Table 1, but with fiN ~ {V2N — 1 + \/2N + l)/2. 



Our interest is primarily with GOE, for which software for exact computa- 
tion appears to be as yet unavailable. Table 3 shows Monte Carlo simulations 
of -P{a;(i) < fj-N + T]\fSia} for quantiles sia of the Fi limit, based on i? = 10^ 
replications. The corresponding 95% confidence intervals have half width 
2^/p^jT^-I>aj X 10"'^ which decreases from 0.001 at pa = 0.5 to 0.0002 at 
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Pa = 0.01 and 0.99. Thus the tabulated values should be correct to within 
±0.001. 



N + 1 


0.01 


0.05 


0.1 


0.3 


0.5 


0.7 


0.9 


0.95 


0.99 


2 


0.010 


0.045 


0.090 


0.279 


0.483 


0.698 


0.914 


0.963 


0.995 


5 


0.012 


0.053 


0.103 


0.300 


0.500 


0.704 


0.907 


0.956 


0.993 


10 


0.011 


0.053 


0.103 


0.302 


0.502 


0.703 


0.904 


0.954 


0.992 


25 


0.011 


0.052 


0.103 


0.302 


0.502 


0.702 


0.902 


0.952 


0.991 


50 


0.011 


0.052 


0.102 


0.301 


0.501 


0.701 


0.901 


0.951 


0.991 


75 


0.011 


0.051 


0.102 


0.302 


0.501 


0.701 


0.901 


0.951 


0.990 


100 


0.010 


0.051 


0.101 


0.301 


0.501 


0.701 


0.901 


0.951 


0.990 


200 


0.011 


0.051 


0.101 


0.301 


0.501 


0.701 


0.901 


0.951 


0.990 


500 


0.010 


0.050 


0.100 


0.300 


0.500 


0.700 


0.901 


0.951 


0.990 



Table 3 

GOE approximation. Let quantiles sia be defined by F-i_{sia) = a for values of a shown 
m the top row, and /in and tn be given by (1). Based on R — 10® replications drawn 
from GOE, table entries are the fraction of replications o/s(i) = — ^n)/tn 

satisfying S(i) < s\a- 



Two features of the numerical results deserve note. First the approxima- 
tions are somewhat better in the near right tail than in the left. This is 
presumably because the underlying approximation of Hermite polynomials 
by the Airy function is anchored at the turning point of the Airy equation 
A"{s) = sA(s), which lies in the right tail at about the 83rd percentile of 
Fi and the 97th percentile of F2. 

Second, the errors in Tables 1 - 3 all have the same sign, suggesting that a 
further shift in the approximating distribution might improve accuracy. We 
experimented in GOE with small changes of the form, setting A^-|_ = A^+1/2, 

(3) /i^(7) = (2iV+ - lN-''^f'\ tm{c) = 2-i/2(iV + c)-l/^ 

and obtained good results. Table 4 and Figure 1, for 7 = 1/5 and c = 1. 

These values differ from //jv and r^v of Theorem 2 by by relative errors of 
0{N~^^^) and 0{N~^) respectively and so have no effect on the validity of 
Theorem 2. However, they provide a substantial numerical improvement, es- 
pecially in the right tail for values of N below 10. Indeed, for some purposes, 
the approximation in the right tail would be adequate even for N = 2. 

Outline of proof. We use the operator norm convergence framework de- 
veloped in Tracy & Widom (2005); our focus of course is on achieving the 
second order convergence rate results. We use the Fredholm determinant 
representations for the finite and limiting distribution functions in terms 
of the two point correlation kernels. A bound of Seiler-Simon, along with 
its orthogonal case analog, bounds the difference in Fredholm determinants 
in terms of the kernels. In turn, the kernels have integral representations 
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TW percentiles 



Fig 1. Left panel: estimated density function for s'f^j^-^ = (^(i) ~ A*jv(l/5))/''"jv(l), compare 
(3), based on R — 10^ samples from GOE2, compared to the Tracy-Widom density function 
fi. Right panel: probability plot of percentiles of the Fi distribution on horizontal axis 
versus ordered values 0/ S(i) on the vertical axis based on the same samples from G0E2. 
Vertical lines mark the 1st, 95th and 99th percentiles of F\ . 



N + 1 0.01 0.05 0.1 0.3 0.5 0.7 0.9 0.95 0.99 

2 0.022 0.073 0.127 0.319 0.505 0.696 0.897 0.950 0.991 

3 0.018 0.067 0.120 0.315 0.505 0.699 0.899 0.951 0.991 

4 0.017 0.063 0.116 0.311 0.505 0.699 0.900 0.951 0.991 

5 0.015 0.061 0.114 0.310 0.504 0.700 0.901 0.951 0.991 
10 0.013 0.056 0.107 0.305 0.502 0.700 0.901 0.951 0.991 
25 0.011 0.053 0.104 0.302 0.501 0.700 0.901 0.951 0.990 
50 0.011 0.052 0.102 0.301 0.500 0.699 0.900 0.950 0.990 
75 0.011 0.052 0.102 0.302 0.500 0.700 0.900 0.950 0.990 
100 0.010 0.051 0.101 0.301 0.500 0.700 0.900 0.951 0.990 
200 0.011 0.051 0.101 0.301 0.500 0.700 0.900 0.950 0.990 
500 0.010 0.050 0.100 0.300 0.500 0.700 0.901 0.951 0.990 

Table 4 

GOE approximation. As for Table 3, but using s'^j^j defined using /ijv(7) and tm(c) given 

by (3), with 7 = 1/5 and c — 1. 



in terms of weighted Hermite polynomials, and so we transfer bounds on 
convergence of Hermite polynomials to the Airy function to bounds on the 
kernels and hence to bounds on the probabilities. 

Convenient uniform bounds on the convergence of weighted Hermite poly- 
nomials to the Airy function come from Liouville- Green theory, which ana- 
lyzes convergence of the solutions of the second-order differential equation 
satisfied by the Hermite polynomials to those of the equation for the Airy 
function. 

The correlation kernels for finite involve polynomials of both degrees 
A^ and A^ — 1, each with its own Liouville-Green centering njv and tiAr_i. 
The overall centering fi^ for the kernel and distribution function must com- 
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bine un and uat-i appropriately to ensure that the generic 0{N~^^^) error 
terms cancel to uncover 0{N~'^^^) convergence. In the unitary case, simple 
averaging suffices: fj^jy = {un + un-i)/'2. For the orthogonal setting, we use 
a formula expressing the GOEat+i kernel in terms of the GUEat kernel plus 
a rank one kernel, and obtain cancellation of 0{N~^/^) errors from these 
two components. 

The Hermite polynomial approximation results are summarized in Sec- 
tion 2. The unitary proof, in Section 3, is a necessary preparation for the 
orthogonal case in Section 4. 

Reproducible code: Matlab files to produce the figures and tables are 
available at the second author's website. 

Related work. Convergence rate results at 0{N~'^/^) were obtained 
by El Karoui (2006) for LITE, Johnstone (2008) for JUE and JOE, and Ma 
(2011) for LOE. The study of Edgeworth type expansions for GUE and LUE 
was initiated by Choup (2006, 2008), who noted that N^^^'^ terms in these 
expansions can be removed by specific choices of the centering constant. 

The Tracy- Widom limit laws for the largest eigenvalue hold much more 
generally — such universality results are an active subject of research. For 
Hermitian Wigner matrices, see Tao & Vu (2010) and references therein, 
and for covariance matrices Soshnikov (2002) and Peche (2009). 

2. Hermite polynomial asymptotics. The Hermite polynomials, Hk{x) 
in notation of Szego (1967, Ch. 4), are orthogonal with respect to the weight 

2 

function w{x) = on (—00,00). The "oscillator wave functions" are nor- 
malized, weighted versions 

with hk = /i7|(x)e-^'dx = ^2''kl 

Classical Plancherel-Rotach asymptotics for H]\[{x) near the largest zero, 
Szego (1967, p. 201) and Anderson et al. (2010, §3.7.2), establish that, for 

ruN = V2N and tn = 2~^/'^N~^/^, 

(4) {2Nf'^TN^N{niN + stn) ^ A{s), 

where throughout we use A to denote the Airy function Ai. 

We will need to explicitly bound the error in the convergence in (4). 
There is now a substantial literature on asymptotic approximations to Her- 
mite polynomials, using for example the steepest descent method for inte- 
grals (e.g. Shi (2008)), the nonlinear steepest descent method for Riemann- 
Hilbert problems (e.g. Wong & Zhang (2007)) and recurrence relations (e.g. 
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Wang & Wong (2011)). Much of this recent attention has focused on expan- 
sions for Hn{\/2N + 1^) and (j)M{-\/2N + 1^) that are valid uniformly for 
large regions of ^. 

For this work, however, we need more detailed information for ^ = l+a^s 
near 1, and specifically uniform bounds for the error of Airy approximation 
for both (j)N and its derivative that have exponential decay in the variable 
s and rate N~'^^^, cf. Proposition 1 below. We have not found this extra 
detail explicitly in the literature, and since the Liouville-Green discussion 
of Olver (1974, Ch. 11) comes with ready-made bounds for approximation 
error for both (j)]\f and we use this as a starting point for extracting, in 
the Appendix, the specific bounds we need. In this section, we explain just 
enough of the approach to describe the bounds we need. 

The Liouville-Green (LG) approach relies on the fact that Hermite poly- 
nomials and hence (pN satisfy a second order differential equation, 

(5) <p'l,{x) = {x^-{2N + l)}cl)j^{x). 

Rescaling the x axis via x = \/2N + 1^, and setting wiyiO = 4>n[x)^ the 
equation takes the form 

(6) w'^{i) = Klf{OwN{0, 

with 

KN = 2N + l, f{^) = e-l. 

The turning points of the differential equation are the zeros of /, namely 
^± = ±1, so named because each separates an interval in which the solution 
is of exponential type from one in which the solution oscillates. The LG 
transformation introduces new independent and dependent variables C and 
W via the equations 



(7) c^^j =/«). w = y^^ 

More precisely, we take 

(8) {mc"\i) = f^^^)di'. 

The transform W approximately satisfies the Airy equation W"{C,) = i^\[C,^{C)-> 
which has linearly independent solutions in terms of Airy functions, tradi- 
tionally denoted by Ai{K?/^C,) and ^\{k?/^C,). Our interest lies in approxi- 
mating the recessive solution K\{k?/^C,). 
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As described in more detail in the Appendix, the error in the Liouville- 
Green approximation can be bounded and one arrives at 

(9) 4>n{x) = {2Nf/^TN^N{x) = eNr{i){A{4\) + 0{N~^)}. 

Here r(^) = [ClO/Cll)]"^^^ is approximately 1 for ^ near 1, and = 
1 + 0{N~'^). This is, then, a version of (4) with an error term of order 
0{N^^), but with Airy function argument k^^C rather than s. 

We focus on x near ujy = \/2N + 1, that is on ^ near the upper turning 
point ^+ = 1. Introduce the rescaled variable s through ^ = 1 + a^s. To 
more closely match the result (4), we want ctn to be chosen so that the Airy 
function argument 

(10) 4^\(1 + aNs) « s 

for s in a suitably large range. A Taylor expansion of the left side yields 

(11) 4\C{l) + aNC{l)s+lals'C)- 

Since C(l) = and C(l) = 2^^^, as follows from (7) and (8), we obtain (10) 
by any choice of the form 

For such a choice, (76) shows that 

+ aNs) = s + 0{s^N-^/^). 

Thus, to replace A{kj^ () in (9) by A{s) entails, in general, accepting an 
error term of 0(iV~2/3) instead of 0(A^~^), and so we use this error scale 
henceforth. 

With the specific choice ctat = tm/um^ we will show that for s > sl, 

l^NiuN + STn) - A{s)\ < CiV-2/3g-s/2_ 

Figure 2 shows that, for values of s corresponding to the bulk of the support 
of Fi, the approximation is tolerably good even for N = 2. 

In fact, since the two point correlation functions depend on both (pN and 
4>N-i-, we need such approximations both for and for (j)N-i- For (pN-i = 
{2N)^/^tn<I)n-i ) the corresponding turning point is at un-i = \/2N — 1, 
though we still use the same scale factor ttv- We use the notation (pNj, with 
Nj = N OT N — 1 respectively, to refer to both cases. In addition, for the 
GOE case, bounds on the convergence of the derivative is also required. In 
the Appendix, we establish 
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Fig 2. Plots of <j>N{uM + tms) for TV = 2, 20 compared with the Airy function A{s) 

Proposition 1. Let sl G K. For s > sl, we have 

(12) \^Nj{uNj+STN)\<Ce-', 

(13) \^Nj{uNj + stn) - A{s)\ < CN~^/\-'/^, 

where the error bounds are uniform in s > sl and N > Nq{sl)- The same 
bounds hold, with modified constants C, when (pNj and A are replaced by 
TNfpMj o-nd A'. 

We record here also some corresponding exponential decay bounds for the 
Airy function and its derivatives A^^\ Indeed, given sl, there exist constants 
Ci{sL) such that 

(14) \s''A('\s)\<a{sL)e-^ s>SL, i = 0,l,2. 

Proposition 1 provides good Airy approximations for both cp^ and 4'N-i, 
but with differing centering values, un and un-i respecitvely. To obtain 
scaling limits for the correlation kernels in GUE and GOE, we need to 
combine these centerings in a manner appropriate to each case. 

It is convenient to express these centering shifts in the rescaled variable 
s. Thus, set 

^^^^ (t)T{s]k) = 4)n{un + tn{s + k^N)), 

tpr{s; I) = (j)N-l{uN-l + mis + /Aat)), 

where 

A;v = [uN - un-i)/tn = N-^'\l + 2-5iV-2 ^ o(iV-^)), 
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— indeed 2^/^A2 = 1.0080! We obtain extensions of Proposition 1: indeed 
from (12), 

\(l)r{s;k)\ < Ce-(*+'=^^) < Ce'', 
and a similar bound holds for |V't('S; k)\. From (15) and Proposition 1, 

(j)r{s; k) = A{s + kAN) + 0{N-^/^e-'^^) 

(16) = A{s) + kAjvA'is) + 0(iV-2/3e-^/2), 

since i(A;Ajv)2|^"(s*)| < CiV-^/Sg-^ using (14) and \s* - s\ < kAj^. 
More generally, but by identical arguments, for r = 0, 1 we have 

(17) #(s; k) = A^'-^s) + fcA7vA('-+i)(s) + 0{N~^/^e-'/^), 

and correspondingly for ijji^\s;k). As a byproduct of a steepest descent 
analysis for Laguerre polynomials, Choup (2006, 2008) derived a three term 
asymptotic expansion for (pj^ and 4>n~i whose first two terms agree with 
(16), though without the uniform error bounds in N and s of Proposition 1. 

3. Unitary Case. Proof of Theorem 1. The argument has three 
components: first we recall determinantal representations of the eigenvalue 
probabilities Fi^^2ixo) and limiting value i*2(so)) along with integral repre- 
sentations of the associated correlation kernels. Second we set up the rescal- 
ing that connects xq and sq, and finally establish the convergence bounds. 

The two point correlation kernel for GUEjv 

7V-1 

SN,2{x,y) = ^ (t>k{x)(l)k{y) 

k=0 

has a useful integral representation (Tracy & Widom 1996, equation(57)). 
Sefi 

(18) <t){x) = {2Nf/^(t)N{x), ^{x) = {2Nfl^(l)N-i{x), 
then 

j'OO 

(19) SN,2{x,y) = \ [(t>{x + z)i){y + z) + i,{x + z)(l){y + z)]dz. 

Jo 

The distribution of may be expressed as a Fredholm determinant 

(20) Fn,2{xo) = P{ max Xfc < xq} = det(/ - S'Ar,2Xo)- 

I'^k'^N 

'^note: our definitions differ by a factor \/2 from those of Tracy and Widom 
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where Xo{x) = I(xo.oo){x) and the operator Sj\f,2Xo is defined via 

/•CO 

{Sn,2Xo)9{x) = i SN,2{x,y)g{y)dy. 

Equivalently, we may speak of 5Ar^2 as an operator on ^2(2:0, 00) with kernel 
SN,2{x,y). On this understanding, we drop further exphcit reference to xo- 
Now change variables, setting x = t{s) = + tns, with //at yet to be 
determined, and xq = r(so). Set also 

(21) Sr{s, t) = TnSn,2{IJ-N + TNS, ^IN + TNt). 

Defining 

(22) (t)r{s) = TN(j){fJ-N + Tjys), Vr(s) = TN^I^il^N + Tns), 

it is clear that (19) becomes 

poo 

(23) Sr{s,t) = \ [(l)r[s + Z)'4jr{t + z)+'4)r{s + Z)(l)r{t + Z)]dz. 

Jo 

Since Sn,2 and Sr have the same eigenvalues, det(/ — S]\f,2) = det(/ — Sr), 
and so 

(24) P{{ma.xxk - iin)/tn < so} = det(/ - Sr). 

Tracy & Widom (1994) showed that the limiting distribution F2 also has 
a determinantal representation 

F2(so) =det(/-5A), 

where Sa denotes the Airy operator on L^(so,oo) with kernel having the 
form 

;>oo 

(25) SA{s,t)= A{s + z)A{t + z)dz. 

Jo 

To derive bounds on the convergence of -F/v,2(a^o) to F2{so), we use a bound 
due to Seller & Simon (1975): 

(26) |det(/-50-det(/-S^)| < - exp(||S,||i + + 1). 

Here ||-||i denotes trace class norm on operators on L^(so, 00). This bound re- 
duces the convergence question to study of convergence of the kernel Sr{s, t) 
to Sa{s, t). 
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Given functions a and b, denote by a o 6 the operator having kernel 

POD 

(ao6)(s,t)= / a{s + z)b{t + z)dz. 
Jo 

In this notation, the kernel difference becomes 

Sr — Sa = (0r O V'r + V'r O f^r) — A O A. 

To facilitate convergence arguments, we rewrite this as 

8(5^ - Sa) = {<Pr +ijT + 2A) O {(l)r + ^Ar - 2A) 

+ {(t)r +A- 2A) O {^r +A + 2A) - " V'r) <> {(pr " V'r)- 

Recall that the centering constant fiN was left unspecified in the defini- 
tions of (pr and ipr in (22). We now choose fi^ so that each term in the preced- 
ing decomposition is 0{N~'^^^). This amounts to choosing the shifts k and 
/ in (15) to satisfy two constraints. First, the centerings fijy = un + kr^^N 
and iJ-N = U]\f^i + /ttvAtv must agree, so that necessarily I = k + 1. Second, 
the N~^/^ term must drop out in the expansion for (p^- + given by (17), 
so that I = —k. We therefore must have, for the present unitary case, 

^27) Ms) = Ms; -I) 

ipris) = Ms; 5) 

which entails that /i^r = (un + ujv„i)/2 as was used in Table 2. Prom 
Proposition 1 and the succeeding discussion we obtain 

Corollary 1 (Complex Case). Let (pj- and ipr be defined by (27) and 
(15). Given sl G M, there exists C = C{sl) such that for N > N{sl) and 

S > SL, 

(28) |,^,(s)| , |^,(s)| < Ce-^ 

(29) \Ms) - A{s)\ , \Ms) - A{s)\ < CN~^l^e-'/\ 

(30) \(Pr{s) + Ms) - 2A{s)\ < CiV-^/Sg-^/a, 

We will need some simple bounds for certain norms of aob. In the unitary 
case, we need the trace norm of a o 6 as an operator on L^(so,oo). In the 
orthogonal case, we need the weighted L^-spaces L'^ {{sq, 00), p{s)ds) and 
L^((so, 00), /9~^(s)(is) for a weight function p such that the reciprocal G 
L^(M). Purther details are given in Section 4. Por some 7 > 0, let 

(31) />(.) = e^N. 

In this section 7 = 0, while values of 7 > will be specified later for GOE. 
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Proposition 2. Let weight functions pi,p2 be chosen from {p,l/p}, 
where p is given by (31), and consider the Hilhert- Schmidt norm of operator 
aob : L?'[p2) — > L^{pi). Assume that, for s > sq, 

\a{s)\ < aNC-^^', \b{s)\ < bNe'''''. 

-(/O < 7 < 2 min(ai, then 

(32) \\aob\\HS < c^^e-('^i+''^)^"±^l^°l, 

ai + Oi 

where C = C{ai,bi,j) = [(oi - j/2){bi - 7/2)]-V2. jf ^ ^f^^^ ^j^^ 
trace norm satisfies the same bound. 

This is a special case of Johnstone (2008, Lemma 7). In the present unitary 
case, we apply Proposition 2, with 7 = 0, to bound the trace norm of each 
term on the right side, using (28) - (30). For each of the three terms, we find 
that ai + 61 > 1 and aNbN < CN""^/^, so that 

ll^r - SaWi < CN~^l^e~"-K 

We may similarly conclude that 

||5,||i <Ce"2^o, IISaIIi < 

Indeed, bounds (28) and (14) show that in each case, and 67V ^ C{sq), 
and that ai = 61 = 1. 

Combining the two previous displays with Seiler-Simon bound (26), we 
obtain Theorem 1. 

4. Orthogonal Case. To establish Theorem 2, we again follow the 
outline of proof given in Section 1. 

1°. Assume that iV + 1 is even. Tracy & Widom (1998) gave a derivation^ 
of the determinant representation 

(33) P{^^max^^Xfc < xq} = A/det(/ - Kn+iXq)- 

Here i^Ar+i is a 2 x 2-matrix valued operator 

(34) KN+i{x,y) = {LSN+i,i){x,y) + K'{x,y), 

^Sinclair (2009) extended Tracy and Widom's derivation to cover + 1 odd, but we 
do not pursue this here. See also Forrester & Mays (2009). 
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where 

(3, -(-)-U-, 2)- 

Here 82 denotes the operator of partial differentiation with respect to the 
second variable, and ei the operator of convolution in the first variable with 
the function e(x) = ^sgn(x). Thus {eiS){x,y) = f e{x — u)S{u,y)du. Finally 
T denotes transposition of variables TS{x,y) = S{y,x). The scalar kernel 



Sn+1,i{x, y) = ^ (l)nix)(j)n{y) + \l ^ (l)N{x)e(t)N+l{v), 

n=0 

and Adler et al. (2000) observe that it may be rewritten as 

(36) Sn+i,i{x, y) = Sn,2{x, y) + 5(/)(x)e^(y), 

where cj) and ij) are as defined at (18). The orthogonal kernel is thus expressed 
in terms of the unitary kernel and a rank one remainder term. The formula 
allows convergence results from the unitary case to be reused, with relatively 
minor modification. 

2°. The limiting distribution has a corresponding determinantal repre- 
sentation 

Fi(so) = ^Jdet{I -Kgoe). 

To state the Tracy & Widom (2005) form for Kgoe, and for the convergence 
argument to follow, it is helpful to rewrite expressions involving e in terms 
of the right-tail integration operator (eg){s) = g{u)du and for kernels 
A{s,t) in the form {eiA){s,t) = A{u,t)du. This is due to the oscillatory 
behavior of the Airy function in the left tail. We write A<^B for the operator 
whose kernel is A{s)B{t). The Tracy-Widom expression states that 

(37) Kgoe{s, t) = (^^^^^ *^|^ _ ^sl^^'^^) 

and the entries of Kgoe are given by 

S{s,t) = {SA-\A®~eA){s,t) + \A{s) 

(38) SD{s,t) = -d2{SA-\A(^~eA){s,t), 

IS{s,t) = -h{SA - \A®~eA){s,t) - \{eA){s) + \{lA){t), 

where Sa is the Airy kernel defined at (25). 
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Defining operator matrices 
we may rewrite (37) in the form 

(39) Kgoe = L{Sa -\A® eA) + \LiA{s) + \L2A{t) + K^ 

3°. We turn to a linear rescaling of formulas (33) and (34). We again set 
X = t{s) = ^n+tns and y = r{t) = HN + TNt, but now with fij\r = fi^ to be 
determined anew in this orthogonal case, see 4° below. Define (pr and ipr as 
before by (22) we occasionally write (f>^ and tp^ to emphasize the different 
centering. We have 

S^{s,t) :=TNSN+l,l{T{s),T{t)) 

= TNSN,2{T{s),T{t)) + irjv0(T(s))(eV')(r(t)) 

(40) =Sr{s,t) + ^Ms){ei^rm, 

where we have used (e-0)(y) = {eipT){t) for a linear rescaling. 
Now det(/ — i^Tv+iXo) = det(/ — Kr)-, where 

Kr{s,t) = TNKN+liT{s),T{t)) 

= TM SN+l,l{r{s),T{t)) + TMK\T{s),T{t)) 



where K"^ was defined at (35). Since det(/ — Kr) is unchanged if the lower 
left entry is divided by tn and the upper right entry multiplied by r^v, 

(41) det(/ - Kn+iXo) = det(/ - Kr) 
where Kr is an operator with matrix kernel 

(42) Kr{s,t) = {LS^){s,t) + K'{s,t). 
Now we rewrite using e and ei. First, define 

/"OO /"OO 

(43) /3n-i = 11 ^r = \{2N)^'^ I <Pn-i, 

-OO 
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and observe that eipr = f^N+i — ^ipr- Thus 

LS^ = L{Sr - \(t>r ® eV'r) + \PN~lL{<i)r ® !)• 

Now L = L + i ^ ^ ^1 and 2(ei + ei ) amounts to integration over M in 
Vei + ei Oy V i i; & 

the first slot. From (23), after interchanging orders of integration and using 

J (j)T = 0, we obtain 

oo POD 

Sr{s,t)ds= / PN-lfl^rit + Z)dz = /3N-ie(f)r{t), 
-oo Jo 

and then 

[(L - L){Sr -4>r(^ eV'r)]2,l = ® e4>r 

as the only non-zero entry of the matrix on the left side. Combining the last 
two displays with (42), we get 

(44) Kr = L{Sr - ® e^) + ^/3n-i[LiMs) + L2Mt)] + K\ 

4°. We now look at the (1, 1) terms in (39) and (44) in order to see, some- 
what informally, how the choice /i^ = un leads to 0(A^~^/^) convergence. 
Thus, we examine the difference 

(45) [Sr - \(t)r ® Cr] " [S A - \A ® lA] + i[/37V„i0^ - A]. 

From definitions (15) and expansions (17), this choice of corresponds to 

c/>f(.) = <A.(s;0) = A(.) + o(iv-2/3), 

^^{s) = Ms- 1) = A{s) + ^N^{s) + 0(iV-2/3). 
We write Ajsi = A + Aat^' and define 

Sa^ = \{AoAn + AnoA). 
From representation (23) and (46), Sr = Saj^ + 0{N~'^^^), while the identity 

Sa^ = Sa- \^nA A 

follows from 



{A<^A' + A' <^A){s,t)= / —[A{s + z)A{t + z)]dz = -A{ 

Jo dz 



s)A{t). 
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Thus 

Sr = SA- A(^A + 0(iV-2/3). 

Since eA' = —A, we have cAn = eA — An A, and so 

^(j)^ ,g) iilj^ = eA - ^An A <S) A + 0{N-'^/'^). 

Forming the difference of the last two displays, we see an important can- 
cellation of the 0{N~^/^) terms involving Ajv, and hence that the first two 
terms of (45) together are 0{N~'^/^). 

A computation with the recursion relation for Hermite polynomials and 
then Stirling's formula (with its 0{N~^) error term) shows that, as — t- oo, 

^^-^"l— J 2(^-i)/2((Ar-l)/2)!-^ + ^^^ 

From this and (46), it follows that the final term of (45) is also 0{N~'^/'^). 

5°. To prepare for the convergence argument for the 2x2 matrix kernels, 
we combine (39) and (44). Noting also from our considerations above that 

SA-\A®eA = Sai, -\A® IAn. 

we obtain the basic difference representation 

- Kgoe = L{Sr -Sa^)- \L{(l)r ^e^jj^-A® IAn) 

+ \Li[PN-l<lyr- A]{s) + \L2[PN-l^r-A]{t), 

from which we may expect to show 0{N~'^/^) convergence, in view of the 
fact that 4)t-,iIJt-,St and I5n-i merge respectively with A^Aj^^Sa;^ and 1 at 
rates of at least 0(iV-2/3). 

6°. We now turn to study the convergence of 
(48) F^+i,i(so) = P{{x(i) - 11n)/tn < so} = - K^) 



to Fi(so) = \/det(/ — Kgoe)- Tracy & Widom (2005) describe with some 
care the nature of the operator convergence of i^Ar+i to Kgoe for the Gaus- 
sian finite N ensemble. We adopt their framework of weighted spaces 
and regularized 2— determinants. Thus, let p be a weight function such that 
G L^(R) and all <?;>Ar G L^{p). Write L'^{p) and L^(p~^) for the spaces 
L'^ {{sq, oo) , p{s)ds) and L'^ {{sq, oo), p~^{s)ds) respectively. 

We consider K^ and Kgoe as members of the collection B of 2x2 Hilbert- 
Schmidt operator matrices B = {Bij, i,j = 1,2) on L'^{p) © L'^{p~^) whose 
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diagonal entries are trace class. Note that e : L'^{p) — L'^{p~^) as a con- 
sequence of the assumption that € L^. The specific p that we use is 
defined in (31) with 7 > 0. 

To analyze the convergence oi p^^i = F/v+i,i(so) to poo = Fi{so), we 
note that their difference is bounded by \p%_^_i — P^l/Poo so that we are led 
to the difference of determinants 

(49) \Fn+i,i{so) - F{so)\ < C{so)\ det(/ - K^) - det(/ - Kgoe)\. 

A Seiler-Simon type bound on the matrix operator determinant for oper- 
ators in B is established in Johnstone (2008). 

Proposition 3. For B,B' g B, we have 

I det(/ -B)- det(/ - B')\ < C{B, B')A{B - B') 

where 

2 

AiB) =Y,\\Bn\\l +Y,\\Bij\\2. 
i=l i^ij 

The coefficient has the form C{B, B') = Y^'^=iCij{ticB,tvB')c2j{B, B'), where 
cij and C2j are continuous functions, the latter with respect to the strong 
(Hubert- Schmidt norm) topology. 

Insert the conclusion of Proposition 3 into (49) to obtain 

(50) \Fn+i,i{so) - Fi(so)| < C{so,Kr,KGOE)A{Kr - Kgoe) 
We exploit decomposition (47), which we write in the form 

Kr - Kgoe = 5' + 6^ + 6f + 5f 

to distinguish a term involving integral kernels, 6^ = L{Sr — Sa^^) from 
terms involving finite rank operators. We establish trace norm bounds for the 
diagonal elements and Hilbert-Schmidt bounds for the off-diagonal entries. 
The distinction between the two norms is moot for the finite rank terms 5f , 
so the trace bounds are actually also needed only for the 6^ term. 

For each term, we show < CN~'^/'^, so that A.{Kr — Kgoe) is 

bounded above by CN'"^/^. We have both \\Kr — KqoeWi and ii — 
tiKGOE converging to at 0(iV~^/^) rate, so that C{Kr, Kgoe) remains 
bounded as — )• 00. 

7°. To bound each term 5ij, we need orthogonal case analogs of the 
uniform bounds of Corollary 1, but now for and their integrals and 

derivatives. From Proposition 1 and the succeeding discussion, we obtain 
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Corollary 2 (Real Case). Let (pr and ip^ be defined by (46) and (15). 
Given sl G there exists C = C{sl) such that for N > N{sl) and s > sl, 



(51) \Ms)\<Ce-', 

(52) < Ce-^ 

(53) IMs) - A{s)\ < CN-^/^-'/^, 

(54) iVr(s) - A{s) - AnA'{s)\ < CiV-2/3e-^/2_ 



The same bounds hold, with modified constants C, when (pr,'4'r,A and A' 
are replaced respectively by (j)'^,^'^,A' and A", or when iprjA and A' are 
replaced by eipr^^A and lA' . 

6^ term. For = L[Sr — Sa^]-, we use Proposition 2 to establish the 
needed Hilbert-Schmidt and trace norm bounds for each entry in the 2x2 
matrix. We write 

St - Sam = i(l)r-A)o^^ + Ao {^Pr - An) + {ipr - An) oc^^ + Ano (0^ - A). 

In turn, for d2{Sr — Sam) we replace the second slot arguments tpr, (V'r — 
An), etc. by their derivatives, and for e{Sr — Saj^), we replace the first slot 
arguments {(j)r — A), etc., by their right tail integrals. 

Consider, for example, the first term {pr — A)o'ipT-. We apply Proposition 2 
using (51) and (53) to set 

aN = CN^'^/^, bN = C, ai = i, 6i = 1. 

The argument is entirely parallel when 82 and ei is applied to {(pr — A)o'ipT-, 
and also for each of the second through fourth terms. Thus, if Dij denotes 
any matrix entry in any component of , we obtain 

(55) llAjll < C7Ar-2/3g-3so/2+7l«ol, 

Finite Rank Terms. As Tracy & Widom (2005) note, the norm of a 
rank one kernel u{x)v{y), when regarded as an operator u®v taking L'^(pi) 
to L?'{p2) is given by 

(56) lln^ull = ||u||2,p2ll^ll2,p-i- 

Here norm can be trace, Hilbert-Schmidt or operator norm, since all agree 
for a rank one operator. 

The finite rank terms include ones of the form L{a® eb). We use (56) to 
establish entrywise bounds 
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(\\a0eb\\ \\a0b\\\ (A+B^ A+B^ 
^ ^ Vl|ea®a|| \\eb(^a\\J - \A_B_ A+B. 

where 

A+ = \\a\\ + , B+ = \\b\\ + 
A_ = ||ea||_, B^ = \\eb\\^. 

Indeed, for the entry apply (56) to aij^bij : L'^{pj) — )■ L'^{pi), where 

pi = p and p2 = p~^ ■ On the right, and henceforth, we abbreviate the Lp' 
norms on L'^{p) and Lp'{p~^) by || • ||+ and || • ||_ respectively. 
Let us indicate how this applies to 

-26q = L[4>r ® e(^^ - An) + {(t)r - A)® eA^] 

Consider the first term on the right side — the second term is similar — and 
apply (57) with a = (j)r-,b = ipr — -^N ■ From Corollary 2 we have 

,/.2p<C7(7)e-2^«+^l^°l, 

-0 

Bl = \\e{ilJr - AN)\t < C(7)iV-^/=^e-*°+^l^ol 
and with similar bounds respectively for and B^. Hence 
(58) A±B± < C(7)iV-2/3e-3^o/2+Tl«ol_ 

Turning to the the 6^ ,62 terms, we have 

^.F ^ ( {uNi -A)® I 0\ , F^t ^ (0 1 {um - eA) 

' \-{uN2-eA)(E)l Oj' l^{uNi-A) 

with uni = I3n-i4't and = I^n-i^4>t- Using (57), we find that the norms 
of the terms in the first column of are bounded by ||ti7vi — j4||+ ||1||- and 
||^^Ar2 — e^ll- ||1||- while the norms of the second column of ((Jf')* are bounded 
by the same quantities interchanged. 

From the definitions, and with sq > 0, we have ||1||?_ < 7~^e~'^'^" and 

\\UN1-A\\+ < I/37V-1 - l|||0r||+ + ||</'r - ^Ik- 
Note that I/Sat-i — 1| = 0{N~^). Using also the bounds of Corollary 2, 
\\um - A\\+ < {CN-^e-'" + CiV-2/3g-so/2)g7^o/2^ 
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and lluATi - ^|| + ||1||_ < CiV~^/^e~*o/^. The term \\un2 - ^M- is bounded 
analogously. 

We finally assemble the bounds obtained from (55), (58) and the analysis 
of 6f and only track the tail dependence on sq for sq > 0. Thus (50) is 
bounded by 

(7jY-2/3|-g-3so/2+7So _j_ 

where the second term results from 6i and 62 ■ It is clear that 7 = 1 yields 
a bound CN~'^/'^e~'^^^ . 



5. Appendix: Hermite polynomial asymptotics near largest zero. 

Define new independent and dependent variables ( and W via the equations 
(7), which put (6) into the form 

(59) ^ = {^2^ + ^.(0^ 

where the perturbation term '0(^) = (~^/'^{(P /d('^){(^/'^). If the perturbation 
term "0(0 were absent, the equation (PW/d^'^ = k^^W would have linearly 
independent solutions in terms of the Airy functions Ai{K'^/^(^) and 61(^2/'^ (^). 
Our interest is in approximating the recessive solution Ai(fi;^/'^C)) so write the 
relevant solution of (59) as W2{C) = ^K^'^^^O'^viO- Iii terms of the original 
independent and dependent variables w and ^, the solution W2 becomes 

(60) W2{t ^) = t^'\i){A{t^'^C) + e2(^, ^)}. 

Olver (1974) — hereafter abbreviated as [O] — provides in his Theorem 
11.3.1 an explicit bound for ri{C,) and hence €2 and its derivative. To describe 
these error bounds even in the oscillatory region of A[x), [O] introduces a 
positive weight function E(a;) > 1 and positive moduli functions M(x) < 1 
and H{x) such that for all x, 

(61) |^(x)| < M(x)E-^(x), \A'{x)\ < U{x)E'\x) 
[Here, E~^{x) denotes 1/E(x).] In addition, 

(62) A(x) = 2"^/2M(x)E-^(x), x>c = -0.37, 
and the asymptotics as x — )■ 00 are given by 

(63) E(x) ~ \/2ei^'^', M(x) ~ 7r-^/2^-l/^ and N(x) ~ vf-^/V/^ 
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The key bounds of [O, Thm. 11.3.1] then state, for ^ > and /(^ = fiO/t 

(64) 162(^,^)1 < (M/E)(K2/3c)[exp{^V(C)} - 1], 

(65) Id^e^iC, k)| < K2/3Ar-i/i/2(^)(N/E)(K2/3c) 

where Aq = 1.04. For k^/^C > c, (62) shows that the coefficient in (64) is 
just ^/2A{K'^/^C)■ Here V(C) = V[^^oo]{H) is the total variation on [^,oo] of 
the error control function H{^) = — Jq^^^ \v\~^^'^^p{v)dv. From [O, p. 403] 
we have XoV[^^oo){H) < 0.28 and hence 

(66) exp{— V(C)} - 1 < 1/iV. 

K 

Application to Hermite polynomials. In the case of Hermite polynomials, 
transformed as in (6), the points ±oo are irregular singularities and the 
points ^± = ±1 are turning points. We are interested in behaviour near the 
upper turning point which is located near the largest (scaled) zero of 
Hj\f. Using (8), the independent variable C('?) is given in terms of /(^) by 

(67) (2/3)C^/2(0 = ^Cie - 1)'/' - I log(e + {e - 1)'/'), 
for ^ > 1, and by 

(2/3)(-c)^/2(0 = ^[cos-i e - c(i - a'/'], 

for C < 1. The function ({^) is increasing and on (0,oo) [e.g. [O, p. 399]], 
with (^(^) non-negative and bounded. It is easily seen that — >• oo as ^ — )■ oo, 
and more precisely, from (67), that 

(68) (2/3)C^/2(0 = i (e' - log? - ^ - log2) + O(r'), 
from which it follows that 

(69) C(6 ~ (4^/3)^/^ as ^ ^ oo. 

We remark that C = C(l) is easily evaluated using L'Hopital's rule. From 
(67), as ^ ^ 1, we have C^(0 = (?^ - 1)/C(0 ^ 2/C, so that C = 2^/^. In 
addition, we shall need the function 

which is positive on (0, oo) since C(^) is strictly increasing. Both r(^) and 
r{S^) are continuous on [0, oo), and as ^ — )■ oo we have r(^) ~ (2^/3)~^/^ and 
r(^) ~ ci.^"''/^, so that r(^) and r'{^) are both bounded on [0,oo). 
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Bound (64) has a double asymptotic property in ^ and k which will be use- 
ful. First, suppose that and hence k are held fixed. As ^ — t- oo, V(C) 
and so from (64) and its following remarks ^2(^,1^) = o{A(n'^^^Q). Con- 
sequently, as ^ — )■ 00 we have W2{S,,i^) ~ C~^/'^{S,)A{k^/'^Q. If the weighted 
polynomial WNiO is a recessive solution of (6), then it must be proportional 
to W2, so that WNiO = cnW2{C, k)- Now cat may be identified by comparing 
the growth of W]\f{C) as ^ — )• 00 with that of ^2(^5 /*) (Section 6): 

(70) c^ = e^"/^4/6(2/Ar)V4^ 

where 6" = 0(1). Now we can use (60) to write (j)N{x) = wj\[{^) in terms of 
the Airy approximation. Below, we write e^v for any term that is uniformly 
1 + 0(A^~^). Hence 

{2Nf/UN{x) = e2'/\]l''w2{C,K) 

Set N+ = N + l/2 and f^v = 2-i/2^r^V6^ g-^^^^ 2'^/^k](^C{0~^^^ = 
and using the Airy approximation (60) to W2{£,,k), we finally have 

(71) {2N)^/^fN^N{x) = CNrmM^TO + e2(e, «7V)}. 

Approximations at degree N and N —1. The kernel SN,2ix, y) is expressed 
in terms of the two functions (j)]\f~i{x) and 07v(x), which need separate 
Liouville-Green asymptotic approximations. Thus, for example, in compar- 
ing the two cases, we have ktv = 2A^ + 1 and kat-i = 2A^ — 1. The turning 
point = 1 and the transformation (^(^) of (67) are the same in both cases, 
hence so is r{^). The analog of (71) then states 

{2N - 2)^/^fN^ict)N^i{x) = ejv-ir(e){^(4^!iC) + e2(e, «7V-i)}. 

Rather than fj\ij = 2"^/"^ N_j_^^^ , we will use the single factor tn = 
2~i/2jv~^/^ in the work below. Clearly, we may replace both (2A^)^/^r7v 
in (71) and (2A^ — 2)^/*f]sf-i in the preceding display by (2A^)-^/^rAr at cost 
of multiplicative error terms e^j = ^ + 0{N~^). 

To summarize then, with the convention that quantities with subscript 
Nj differ for Nj = N, N — 1, while those with subscript A^ do not, we have 

(72) i2N)'/^TNC^Nj{x) = eNjrmAiKjijC) + ^2{i, ^^Nj)], 

Denote the left side of (72) by (pNj- We seek a uniform bound on the Airy 
approximation. If we write x = .Jk^^ in the form u^j + stat, then we have 
in particular un = \/2N + 1 and U]\f-i = \/2N — 1. In turn 

C = 1 + STN/y/K:^j = 1 + SaNj, 
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where we define 

(73) aNj = tn/unj = 2~^/^N~^/%2N ±1)-'^/^ = 2-'^N-^/^{l + 0{N~^)). 

We turn now to the proof of Proposition 1. We first record some properties 
of the map s — )■ k^^C(1 + (^Njs), which we sometimes abbreviate as k^/'^C- 

Lemma 1. Given sl G K, 

(74) 1^2/3^ _ s| < |s|/4 for SL<s< iV^/^ N > Nq 

(75) 1^2/3(^1 < c\s\/4 for sl < s < CN^'^, all N. 

Proof. Expand C(C) about the turning point = 1: 

(76) /^^j C(l + (^NjS) = K^^jCTNjCs + ll<^Nj(^NjS'^CiC*)- 

We note from the definitions that 

(77) 4>iv,C = (1 ± l/{2N)y/^ = 1 + 6^, 

with \5n\ < for all > 1. Since < (" is bounded, we find that 

P«) i-'^'^c -.1^(1 + ^) 1.1. 

again for all > 1. If s < N^^^ then the right side is bounded by |s|/4 for 
N > No{M,sl). If |s| < N^^^ then we have (75) for C = C{sl, M). □ 

We consider some global bounds, valid for s > sl, or equivalently for 

> 1 + SLfJNj- 

Lemma 2. Let sl < 0. Let ^ = l + a^jS with cjMj satisfying (73). There 
exists C = C{sl) such that for s > sl, 

E"H4•C)<Ce-2^ 

N(k2/3(^) < + |s|l/3). 

Some immediate consequences: using (61) and M < 1, for s > sl, 

(79) |^(4-C(0)I < |M/E|(k2/3c) < Ce-2^ 

(80) |e2(C,'^)l < N~^\U/E\{k^/^C) < CN~^e-^' 

(81) \A'{4^Cim < |N/E|(k2/3c) < ai + \s\^/^)e-^^ 
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Proof. First, since /(^) = (^ + 1)(^-1) > 2aNjS, we use (77) to observe 
that for s > r^, 

KNj(^Nj^/f > ^K'NjO-%'^j\fs > eNjr. 
Hence, from (8), again for s > r^, 

Now choose r large enough so that for N > A'^o and j = iV, — 1 , we have 
GNjT ^ 1- From (63) we have E~-'^(s) < Cexp(— |s'^/^) for s > 0, and so in 
particular for s > r^, 

E-^(4;0<C(r)e-^ 

For s G [si,,r^], we simply use the bound E > 1, 

For the second statement, we will use the bound N(s) < 1 + ([O, 
p. 396-7]). First, for s < A'"^/^, using the bound on N and (75), we obtain 
N(k2/3() < C(l + |s|^/^). When s > N^/^, we use (67) to bound 

{2/?,)C'l\i) < j\dt< e/2 < 1 + al^s"" < coaljs\ 

From (73) we have K,NjC^Nj = 2N~^/^ and so unjC^^'^ < cis^ for aU N and 
hence N(k;^/^C) < 1 + c\^^s^/^ as required. □ 

Proof of Proposition 1. We begin from the formula 

(82) (l)Nj{uNj + stn) = CNjriOiMl^Nj C) + e2(C, I^Nj)}- 

The bound (12) then follows from (80), (79) and boundedness of r(^). To 
ease notation, we will as needed drop subscripts from eNj,(T]\fj, knj, and tn, 
writing e for a term that is generically 1 + 0(iV~^). 
For the next bound, we differentiate (82), obtaining 

(83) TNCp'jyj [uNj + stn) = Di + D2 + D3, 

with the component terms given by 

Di = ear(e)[^(K2/3C) + 62(^,^^)1 
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Since r(^) is bounded, we again use (80) and (79) to conclude that < 
C(TAre-2^ for all s > sl- From (77) and (69), we observe that 

(84) cjn4'C{0 < C\e^' < C(l + ali'\s\y-'). 
Turning to the second term, we have from (84) and (81) that 

\D2\ < C\A'{k^/^C)(^k^/^C{0\ < C(l + a]l'^\s\'^/^){l + \s\'^/^)e~^' < Ce-^ 
Using (65) and (66), we can rewrite as 

m < CN~' . r(OfT7v4'C(0 • (N/E)(«:2/3<^). 
Using also boundedness of r(^), (84) and (81), we conclude for all s > sl, 

l^sl < CN-'^ • V2 • (1 + a]l'^\s\'^/^){l + \s\^/^)e-^' < CN'^/'^e-'. 
This completes the proof of bound (12) for T]\f(f)]\fj. 

Turning to the error bound (13) and its analog for T]\j(j)i^j, we first note 
that we may confine attention to s G [sl, N'^^^], since for s > N^/^, the 
bounds follow trivially from (12) and its analog and (14). 

We use the decomposition suggested by (82): 

^Nj{x) - A{s) = [eN,m - m4jO + [m4jO - Ms)] 

= Eni + En2 + En3- 
For the E^i term, first use ^ = 1 + cr^jS to write 

(85) |C(e)/C - 1| = / C{u)/C du < CsGNj. 



since C,{u) is bounded for u G [1, 1 + sa^j] C [1, 1 + ^/^]. Together with 
r{i) = m)K]-^'\ this yields 

(86) \eNAi)-^\<C{l + s)aNj. 

Combined with (79), we obtain |^Ari| < CaNj{l + s)e-^' < CA^~2/3g-s/2^ 
For the En2 term, we use (74) and (78) to write 

\A{4]C)- A{s)\<C\s\N-^l\N-^/^ + \s\)snv{\A'{t)\ : ls<t<\s] 
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uniformly for s G [sl, N^^^], where we used (14). 

Finally, for the -Bats term, (80) and boundedness of r imply that 

\En3\ < e^,-r(e)(M/E)(4/|c)iV~' < CN~^e"^'. 

We turn now to the proof of (13) for TN4>Nj on [sl, N^^^]. Using (83), we 
may write the difference TN4>'j^j{x) — A'{s) as Di{s) + [D2{s) — A' {s)] + D^ls) . 
Now decompose D2{s) — A'{s) = Gi + G2 + G3, with 

Gi = MO - i][m/c]A'{K'/'c), G2 = m/c - ma^k^/'^o 

and G3 = A'{k'^/^C) - A'{s). Combining (86), (81) and then (85), we find 
< C(l + s)aN • 2 • C(l + \s\'^/^)e^^' < CiV-^/Sg-^/^ 
IG2I < CsaN ■ C{1 + \s\'^/^)e-^' < CiV-2/3g-./2^ 

G3 is treated in exactly the same manner as the Ei\f2 term above, additionally 
using the equation A"{x) = xA{x). 

6. Identification of cjv We first remark that as C — 00 when ^ — t- 00, 
we may substitute A{x) ~ [2^yl^x'^/^]-'^ e^p{-{2/3)x'^/'^} into (60), along 
with (-^/^ = [C/fiOV^^ from (7) to obtain 

Consequently we may express cat in terms of the limit 

= hm wn{0 ■ 2^K^/^f^'\0 ■ exp{(2/3)KC=^/2|^ 

Write Nj^ for + 1/2. Since wn{0 = 't'Nix) with x = -y/SiV+C, and since 
H]\f{x) ~ 2^x^ , we have as ^ — )• 00, 

WNiO = h-"\~^^^'H^{^i) ~ /.-^/^e-^+^^+^'°s«2^(2iV+)^/^ 

and f^/*{i) = e(l°g?)/2+0«"'), while from (68) 

exp{(2/3X3/2} = giV+5^-^+log?-7V+/2-iV+log2+0(r^). 

Multiply the last three quantities: the coefficients of and log.^ cancel, 
leaving ^—dependence of only 0(^~^) as ^ — )■ 00. Hence 

CM = 2V^A^i/6/i^i/2(2iV+)^/2e-^+/2 2^-^+, 

and noting that {N/2) log iV+ = {N/2) log iV + 1/4 + 0(iV"^), we get 

K-^'^CN^/h^ = \/2^exp{f log2 + f logiV- f + 0(^)}, 

Applying Stirling's formula to hiq = \/ti2^ N\, and dividing into the previous 
display yields (70). 
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